Ultrafast Molecular Imaging by Laser Induced Electron Diffraction 
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We address the feasibility of imaging geometric and orbital structure of a polyatomic molecule on 
an attosecond time-scale using the laser induced electron diffraction (LIED) technique. We present 
numerical results for the highest molecular orbitals of the CO 2 molecule excited by a near infrared 
few-cycle laser pulse. The molecular geometry (bond-lengths) is determined within 3% of accuracy 
from a diffraction pattern which also reflects the nodal properties of the initial molecular orbital. 
Robustness of the structure determination is discussed with respect to vibrational and rotational 
motions with a complete interpretation of the laser-induced mechanisms. 



With the development of attosecond laser sources 
ultrafast molecular imaging has become a major research 
topic in modern physics. On one hand, attosecond laser 
pulses are directly used to image dynamical processes 
in schemes such as the attosecond pump-probe spectro- 
scopic mapping of molecular motions [2] or the interfer- 
ometric real-time observation of electronic motions 
On the other hand, a number of schemes have been pro- 
posed to image molecular structure, which are all based 
on the rescattering mechanism |4|]. This phenomenon 
consists of the tunnel ionization of an electron followed 
by its acceleration and its return by the field, to end 
with its recollision with the molecular ionic core. The 
possible outcomes of this mechanism are all relevant to 
molecular imaging. The elastic scattering of the return- 
ing electron with the ion core defines the LIED Q, 
which can be compared with ultrafast electron diffrac- 
tion using an external electron source [7] . Inelastic scat- 
tering of the returning electron may be accompanied by 
electronic excitation of the parent ion, its further ion- 
ization or the emission of high-energy radiation. This 
last process, high-harmonic generation (HHG), has been 
used directly as a probe or indirectly to unveil orbital 
and molecular structure, and even to image rotational 
and vibrational motions 0|. 

We wish to assess the feasibility of imaging polyatomic 
molecular structure, e.g., measuring its geometrical pa- 
rameters, such as bond- lengths and/or bond- angles), on 
an attosecond time-scale using the LIED technique. We 
show how a diffraction pattern constructed from mea- 
surable photo-electron momentum distribution encodes 
informations on the electronic orbital nodal properties 
but also on the geometry of the nuclei. We assess the 
robustness of this two-fold structure determination pro- 
cedure with respect to inevitable uncertainties about the 
alignment of the molecule relative to the field polariza- 
tion direction, its bond lengths and angles. We consider a 
relatively simple model of the CO2 molecule, the choice of 
which being motivated by different considerations: Being 
linear in its symmetric equilibrium geometry, it is sim- 
ple enough, in particular with respect to its alignment 
properties, and yet with three nuclei, it has enough in- 
ternal degrees of freedom to make its molecular structure 



determination challenging. 

To concentrate on the problem of reading geometrical 
informations out of the photo-electron diffraction pat- 
tern in momentum space, we develop a theoretical model 
based on an effective single active electron (SAE) Hamil- 
tonian. Nuclear dynamics and its coupling to the elec- 
tron ionized out of a definite orbital (initial state) can 
be ignored at the time scale of a few-cycle pulse. Nu- 
clear degrees of freedom are thus taken into account only 
through statistical distributions of CO bond-lengths R 
and polar angles with respect to the laboratory y-axis 
along which CO2 is supposed initially aligned. The in- 
tense laser pulse, which monitors the electron dynamics, 
is linearly polarized along the laboratory x-axis and the 
essential of the electron diffraction process takes place 
in the (x, y)-plane. Electronic wave packets are gener- 
ated by solving the time-dependent Schrodinger equation 
(atomic units are used): 
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with a SAE Hamiltonian using a soft- Coulomb potential 
and an electron-laser dipole coupling in the length gauge: 
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r and p being the electron position and momentum, and 
a labeling the three nuclei located at fixed position vec- 
tors p a . The laser field S(t) is characterized by the angu- 
lar frequency uol and the amplitude S'q. The r-dependent 
effective charge Z a (r) = + (Z° - Z^°)exp(-|r - 
Pa\ 2 involves a single adjustable parameter <j a . It 
is defined, together with the softening parameter a a , so 
as to give field-free eigenstates of the SAE Hamiltonian 
identifiable with the highest occupied molecular orbitals, 
HOMO— n (n = 0, 1,2), in the sense of giving the same 
ionization energy and the same nodal properties as the 
actual molecular orbital of the 22-electron molecule. Nu- 
merical wave packet propagation uses a third-order split 
operator formula on a finite 2D grid continued analyti- 
cally by projection onto Volkov states for \r\ > 45 A. 

The four panels on the second row of Fig.[T]show the re- 
sults of the time-propagation from an initial state 3>(r, to) 
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FIG. 1. (Color online) Ionized electron motion in laser-driven 
CO2 molecule. First row, from left to right: Field- free (tv 9 ) 
HOMO of CO2, external field's waveform, and calculated 2D 
asymptotic electron momentum distribution for the CO bond 
length 4.8 A. Various snapshots of the calculated 2D elec- 
tron density are shown in panels (a) to (d) in the second 
row, corresponding to different times within the optical cy- 
cle (ujl — 0.06 a. u., ^0 = 0.15 a. u.) as marked on the field 
waveform above. 
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FIG. 2. (Color online) Averaged normalized electron diffrac- 
tion pattern S(k y ) obtained with a single cycle pulse (see 
Fig. [J) for a perfectly aligned molecule, with R — 4.8 A. In 
the upper panel, the black solid and green dashed lines cor- 
respond to the diffraction patterns obtained from the (ir g ) 
HOMO and (a u ) HOMO-2 molecular orbitals of C0 2 . The 
red solid line in the lower panel is for the (tt u ) HOMO-1 initial 
state. 



taken as the HOMO, characterized by a perpendicular 
nodal plane which is preserved through time-evolution, 
when the laser polarization is orthogonal to the molecu- 
lar axis. Each panel depicts the electronic density at a 
different time within a single optical cycle, as indicated 
by the corresponding letter on the waveform shown in the 
middle panel of the first row. The most important ob- 
servations from these spatial wave functions are first the 
fringes of panel (c) obtained during the forward (x < 0) 
electron motion as a signature of the two oxygen atoms 
acting as electron ejection sources, and second the rich 
pattern of interference fringes of panel (d) in the back- 
ward (x > 0) motion resulting from the electron diffrac- 
tion analogous to Young's slit experiment with photons. 

More relevant for structure determination is the 
diffraction pattern in reciprocal momentum space, ac- 
cessible experimentally through time-of-flight electron 
velocity mapping, and determined numerically by the 
Fourier transform <&(k x ,k y ) of the asymptotic electron 
wave packet. To illustrate this, we show, in the right- 
most panel on the first row of Fig.HJ \&(h x , ky)\ 2 f° r the 
CO bond length R = 4.8 A, corresponding to a symmet- 
rically stretched molecule. The features of this momen- 
tum distribution, with a succession of vertical interfer- 
ence fringes, are in agreement with the interpretation of 
Ref. 0. From this distribution, a /^-averaged diffraction 
pattern can be defined as 

S{ky) = f^(k X ,ky)\ 2 dk X , (3) 

where Y denotes the domain of integration. In anal- 
ogy with Young's double slit experiment, we expect a 



good resolution of the interference signal if the electron 
de Broglie wavelength Adb = h/p is smaller than the 
slit separation d = 2R. This criterion corresponds to 
\k\ > tt/R and therefore to the high energy range of the 
photo-electron spectrum. To insure a resolution of the 
order of 0.5 A, we define here the domain of integration 
T as the region where \k\ > 3.15 a. u. 

For a single-cycle pulse, and a perfect alignment sit- 
uation, the averaged distribution S(k y ) associated with 
the diffraction of an electron ejected from the HOMO is 
shown as a black solid line in the upper panel of Fig.[2j 
One can notice that S(k y ) exhibits clear interference 
fringes at large k y , characterized by a regular succession 
of peaks, from which geometrical parameters will be in- 
ferred. Indeed, that S(k y ) directly gives informations 
on the nuclear geometry can be seen by the following 
simple symmetry considerations: The CO2 HOMO ini- 
tial molecular orbital is anti-symmetric with respect to 
the perpendicular ah mirror plane containing the carbon 
atom, and so is the electronic wave function at all time, 
assuming perfect alignment of the molecule orthogonal to 
the laser polarization. It can approximately be written 
as 

(r, t) (x f{x, y-R,t)- f(x, y + R,t) , (4) 

which, through Fourier transformation, gives 

4v 9 (M)oc (e iiifc « -e- iRk y)f{k,t). (5) 

The squared modulus of this wave function behaves as 
sin 2 (Rk y ) and thus exhibits zeros for k y = nir/R. The 
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positions of the dark fringes should thus be directly re- 
lated to the internuclear distance R. These positions are 
indicated by dotted vertical gray lines in Fig.[2j The lo- 
cal minima in the calculated spectrum S(k y ) are seen 
to agree fairly well with these predicted positions. The 
complete quantitative analysis gives, within the model 
assumptions, R ~ 4.92 A, with a deviation of less than 
3% from the initial input value R = 4.8 A. 

Given the role that lower-lying valence molecular or- 
bitals play for CO2 in intense laser fields [9], the averaged 
S(k y ) distribution for the HOMO-2 is shown as a green 
dotted line in the upper panel of Fig.[2j Though a bond- 
ing a u orbital, the HOMO-2 presents the same symmetry 
property as the i\ g HOMO with respect to the ah mirror 
plane containing the carbon atom. Eqs. (|4]) and (|5j) are 
thus valid for the HOMO-2 also, and the photo-electron 
signal exhibits the same interference pattern for large k y 
values. In contrast, the bonding tt u HOMO-1, being sym- 
metric with respect to the ah mirror plane, gives rise to 
fringes out-of-phase with respect to those associated with 
the HOMO and HOMO-2. The nodal structure of this 
orbital is reflected through 

(r, t)<xf{x,y- R, t) + f(x, y, t) + f{x,y + R, t) , (6) 

which, through Fourier transformation, gives 
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The squared modulus of this wave function, which be- 
haves as [1 + 2 cos(Rky)] 2 , is characterized by a regular 
succession of two peaks of different amplitudes. Their 
maxima, located at k y = nir/R, and their relative am- 
plitudes are reproduced in the full quantum calculation 
shown in the lower panel of Fig.[2j The total photo- 
electron signal, constituted by combined contributions 
from the HOMO, HOMO-1 and HOMO-2 which can be 
controlled by the field intensity, could be used to read 
accurately the CO bond length from the fringe pattern 
observed in the large Arrange. 

The structural reading made above of the diffraction 
pattern refers to a perfectly aligned molecule in a fixed 
nuclear geometry configuration. In reality, molecules are 
subject to rotations and vibrations, and these geometri- 
cal parameters are usually distributed over some range 
dictated by the rotational and vibrational state in which 
the parent molecule is initially prepared. A first aligning 
laser pulse can be used to create a rotational wave packet 
which periodically rephases, thus yielding, at the revival 
time, a strong molecular alignment along the polarization 
direction. Angular distributions confined in a cone an- 
gle of about 20° have been reached experimentally with 
CO2 [10]. The left column of Fig. [3] shows, in red dot- 
ted lines, the diffraction pattern incoherently averaged 
over a Gaussian distribution of alignment angles with 
a characteristic width of 20° for R = 1.2 A (first row) and 
R = 4.8 A (second row). This Figure also shows, in black 
solid lines, the result obtained for a perfect alignment. 
The contrast loss seen in this Figure denotes a certain 
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FIG. 3. (Color online) Averaged normalized electron diffrac- 
tion patterns S(k y ) resulting from ionizing the HOMO of CO2 
under the same field conditions as in Fig.filfor R= 1.2 A (first 
row) and R — 4.8 A (second row). The solid black lines cor- 
respond to a perfect alignment at fixed internuclear distance 
for a single optical cycle pulse. The red dashed lines corre- 
spond to, column (a), an imperfect alignment, column (b), a 
distribution of R and, column (c), a 10 fs FWHM sine-square 
pulse. 



sensitivity to alignment defects. However, the fringes are 
still clearly observable at their original positions in the 
range \k y \ > 0.7 a. u. The extraction of geometrical infor- 
mations thus remains achievable even in case of imper- 
fect initial alignment of the molecular target. To study 
the effect of variations of the bond length, we also con- 
sidered a Gaussian distribution of internuclear distances 
with a characteristic width of 0.2 A. Note that this width 
is larger than the typical size of the ground vibrational 
level associated with the stretching and bending modes 
of CO2. The results are shown in the middle column of 
Fig. [3l One can notice that the diffraction pattern S(k y ) 
is almost insensitive to the initial distribution of R. The 
zeroes in the averaged S(k y ) remain basically unshifted 
and can therefore be used for geometrical information 
and imaging purposes. To conclude, we expect that a 
measurement on imperfectly aligned CO2 molecules in 
their ground vibrational level would provide diffraction 
patterns that are no less analyzable than those shown in 
Fig.[2j and from which the CO equilibrium bond length 
could be extracted unambiguously. We also note in Fig. [3] 
that, as expected, the oscillation period of S(k y ) is 4 
times smaller with R = 4.8 A than with R = 1.2 A. 

Finally, the right column of Fig. [3] shows the re- 
sults obtained with a 10 fs FWHM sine-square pulse of 
same carrier wavelength, 800 nm. In this case, the 2D 
photo-electron momentum spectrum (not shown) is dif- 
ferent from the one obtained under the single-cycle pulse 
(Fig.[IJ. It is more symmetric with respect to k x , re- 
flecting a different rescattering physics: Concentrating 
on what happens near the peak of the pulse envelope, one 
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distinguishes two separate electron trajectories, which 
correspond to ionization events taking place at different 
times. The associated electron wave packets are then 
driven by the field in opposite directions. Contrary to 
what happens in the single-cycle case, these trajectories 
can both return to the ion and contribute to the diffrac- 
tion event. The returning electron waves overlap and in- 
terfere with each other in the low-energy region of the 2D 
momentum distribution, corresponding to \k\ < 2a.u., 
blurring out the diffraction pattern and prohibiting a 
simple reading of the fringes in this region. However, 
even with a long pulse duration, only a few trajectories 
contribute to the highest energy domain of the spectrum. 
An average of the 2D diffraction map over the domain T 
of Eq. (|3]) corresponding for instance to Ifc^l > 3.15 a. u. 
helps recovering clearly visible interference fringes. In 
the right column of Fig.[3j the black solid line shows the 
result obtained using this integration domain in the case 
of a single cycle pulse, while the red dotted line refers to 
a 10 fs FWHM sine-square pulse. It is seen that with this 
pulse, not only does the diffraction signal S(k y ) continue 
to reveal the geometrical and orbital structure of the par- 
ent molecule, but it also is amplified by about an order 
of magnitude as compared to the single-cycle signal. 

To conclude, we have suggested a simple yet robust 
method to extract the molecular structure from the 
photo-electron spectra of a laser-driven linear, symmetric 
molecule. It is applicable to any sufficiently symmetric 
molecular system for which the HOMO, or any molecular 
orbital contributing to the ionization signal, has specific 
symmetry properties that would allow for a simple inter- 
pretation of the diffraction patterns. The detailed struc- 
ture of the diffraction image reflects the symmetry of the 
molecular orbital from which the recolliding electron em- 
anates since this symmetry is conserved during the time- 
evolution of the system under the combined effect of the 



Coulomb forces and the laser field. Such a symmetry 
conservation actually goes beyond the SAE approxima- 
tion, and we expect the diffraction images calculated here 
to survive inclusion of many-electron effects, contrary to 
what was found relative to corrections to the strong field 
approximation [14j , which is not assumed here. 

As the acquisition process, relying on the rescattering 
mechanism, is faster than nuclear motions, it can be used 
to produce a sequence of stroboscopic pictures, making 
up, for example, a movie of vibrational dynamics. The 
imaging technique could be implemented using electron 
spectroscopy based on different schemes such as angle- 
resolved time-of-flight detection or velocity map imaging 
pr| . These methods are currently used in the atomic, 
molecular and optical physics community, and offer an 
energy resolution Ae/e of the order or better than 1% at 
50 eV, which turns out to be largely sufficient to resolve 
the fringe separations suggested by the present simula- 
tions. While, in the present theoretical study, a typical 
800 nm laser radiation was considered, preliminary nu- 
merical simulations also show that our proposal could be 
used efficiently with few-cycle pulses produced nowadays 
in the IR domain (1-3 /im) [12j. Practically, a single-cycle 
pulse has been produced at A=1.8 /im [13]. Higher rec- 
ollision energies are then obtained, yielding shorter elec- 
tron de Broglie wavelengths. These laser systems should 
therefore provide a better resolution in terms of internu- 
clear distances at fixed laser intensity. 
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